Search for Correlated High Energy Cosmic Ray Events with CHICOS 
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We present the results of a search for time correlations in high energy cosmic ray data (primary 
E > 10 14 eV) collected by the California High school Cosmic ray ObServatory (CHICOS) array. 
Data from 60 detector sites spread over an area of 400 km 2 were studied for evidence of isolated 
events separated by more than 1 km with coincidence times ranging from 1 //sec up to 1 second. 
The results are consistent with the absence of excess coincidences except for a 2.9a excess observed 
for coincidence times less than 10 (j,sec. We report upper limits for the coincidence probability as a 
function of coincidence time. 

PACS numbers: 96.40.-z, 96.40.Pq, 98.70.Sa 



The spectrum of cosmic rays above 10 15 eV has been 
studied up to energies of 10 20 eV by observing the large 
extensive air showers created by the primary incident 
particles Below 10 15 eV there have been direct mea- 
surements via instruments flown in satellites or as bal- 
loon payloads @ The energy spectrum falls steeply 
with a power law E~ 2J for energies up to the "knee" 
at 4 x 10 15 eV, and even more steeply, E~ 3 , beyond 
the knee. Shock acceleration in supernovae provides a 
successful explanation up to the knee, while the source 
of particles beyond the knee region is speculative. Above 
10 19 - 6 eV, extra-galactic protons should interact with the 
cosmic microwave background (CMB) and thereby lose 
energy, resulting in a sharp decrease in the number of 
cosmic rays with energies above the Greisen-Zatsepin- 
Kuzmin (GZK) limit of 10 19 6 eV 0. Experimental data 
regarding the existence of the GZK cutoff are inconclu- 
sive at this point 0,O|. 

Correlations between cosmic rays would indicate that 
the particles have some common history. Such correla- 
tions could provide information about the source of the 
cosmic rays, the number or distribution of sources, or 
about the propagation of cosmic rays. For example, a 
recent study of ultra-high energy cosmic rays by AGASA 
indicates clustering of the directions of origin in the 
sky. 
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This paper addresses the possibility that isolated cos- 
mic ray events separated by > 1 km arrive in time coin- 
cidence. Such correlated cosmic ray events could result, 
for example, from the photodisintegration of heavy nu- 
clei (i.e., iron) by solar photons Q. A previous study, 
: 8j, indicated episodic evidence for time correlations up 
to 10~ 3 s in events separated by ~ 100 km. A more recent 
search for correlated events j^] at large distance scales, 
~ 500 km, found a few candidate events but these were 
consistent with interpretation as accidental coincidences 
between uncorrelated events. We have studied 17 months 
of data obtained with the CHICOS array during January 
2003 through July 2004 and searched for evidence of cor- 
related air shower events separated by > 1 km with en- 
ergy threshold 10 14 eV. This paper reports the results of 
that search. 



CHICOS 

The California High school Cosmic ray ObServatory 
(CHICOS) observes cosmic ray induced air showers with 
an array of detector sites located on school roofs in the 
Los Angeles area (Lat. 34.1°, Long. —118.1°, average 
245 m above sea level). The sites are separated by dis- 
tances of typically 2 — 3 kilometers, with the overall array 
covering an area of ~ 400 km 2 . As shown in Fig.^ the 
detector sites are located in two major groups in the San 
Gabriel and San Fernando valleys, separated by ~ 35 km. 
These 60 sites contain four pairs separated by less than 
1 km, and all other pairs of sites have larger separations. 
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FIG. 1: Locations of the 60 operating CHICOS sites used 
in the analysis in this paper are shown as open circles. The 
Caltech site is located at the origin, x indicates distance in 
the easterly direction, and y indicates the distance in the 
northerly direction. The sites in the San Gabriel valley are 
clustered around the origin, whereas the sites in the San Fer- 
nando valley are centered about 35 km to the west. The rilled 
circles are not included in the > 1 km analysis. 



During the period corresponding to the dataset reported 
in this paper, the number of operational sites increased 
from 31 to 60. 

Every detector site contains two plastic scintillator de- 
tectors, separated by ~ 3 meters, with each detector 
having ~ 1 m 2 area and 5—10 cm thickness. Photo- 
multipler signals are processed by a custom-built time- 
over-threshold discriminator circuit. Timing of the detec- 
tor signals and the GPS receiver signals is facilitated by 
use of National Instruments 6602 80MHz timer/counter 
cards. We have verified that the GPS receiver, either Mo- 
torola UT+ or M12, provides relative timing accuracy of 
~ 50 nsec. Pulse heights, determined from the time-over- 
threshold measurements, are calibrated and continuously 
monitored using the high flux (~ 200 Hz per detector) of 
incident throughgoing muons. The discriminator thresh- 
old is set at ~ 1/4 of the observed muon peak. 

Data are stored on local hard disk and automatically 
transferred to Caltech via internet every night by the 
computer located at each site. "Trigger" events are de- 
fined as those where both detectors at a site record signals 
greater than 2 single vertical particles. Each site gener- 
ates about 1000 trigger events per day. These "trigger" 
events and all "match" hits (single-detector hits within 
50 /isec of any trigger elsewhere in the array) are trans- 
ferred to Caltech. Shower-building software on our server 
then analyzes these data every morning and constructs 
files of built showers for further analysis. 

Extensive air showers generate coincidences among 
several sites. However, the core of the air shower that will 
generate trigger events is generally much smaller than the 
spacing between CHICOS sites. Therefore, a typical air 
shower event involves one trigger with "match" events 
at the neighboring sites. The trigger events are used 
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FIG. 2: Histogram of the excess of successive trigger pairs vs. 
the time difference. The data show a large signal associated 
with the tails of air showers (solid circles) that disappears 
when the distance between the pairs is required to be greater 
than 1 km (open squares). 



to locate the cores of potential air showers, where the 
particle density is relatively high. Showers that trigger 
one site with matching hits at several neighboring sites 
must have extremely high energies of > 10 18 eV. We do 
observe such large air shower events with the CHICOS 
array, presently at a rate of about one per month. The 
arrival times and intensities of detector signals provide 
the information necessary to reconstruct the incident di- 
rection and energy of the primary particle that created 
the air shower. 

Most trigger events are isolated single events (i.e., no 
nearby matches) which are generated by much more fre- 
quent lower energy showers with a threshold of about 

10 14 eV. The rate of these triggers is comparable to ex- 
pectations based on the previously measured flux p| and 
computer simulations of air showers with the AIRES code 
[10|. In this paper we study double trigger events where 
two sites separated by more than 1 km both record trig- 
ger events within a certain coincidence time. 

In a separate study, we deployed 4 additional sites (not 
shown in Fig. 1) separated by 100 — 400 m on the Cal- 
tech campus to provide us with a more abundant sam- 
ple (~ 10 per day) of lower energy showers in the range 

10 15 — 10 18 eV. These data give us confidence that the 
CHICOS equipment functions as expected, exhibiting the 
correlated data that should result from air showers. In 
particular, the shower properties can be reconstructed 
and the events are well described by single air showers. 
These data are under analysis and will be reported in a 
future publication. We note that in this dataset, double 
triggers are not uncommon features of larger air show- 
ers with a core size capable of generating triggers at two 
sites separated by < 600 m. The rate of these double 
triggers falls steeply with distance between sites, and is 
well described by AIRES simulations. In this paper we 
are studying the possibility of double triggers at sites sep- 
arated by > 1 km, so data from these additional 4 sites 
on the Caltech campus are not used in the analysis in 
this paper. 
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Correlation Analysis 

The main focus of this analysis has been the trigger 
data sequence, which forms a complete record of all the 
events detected by the array with sufficient energy to 
trigger a single site. In order to examine these data for 
time correlations, a randomized data set was constructed 
directly from the real data, as in @ , by offsetting the se- 
quence of triggers at each site by some integer number of 
seconds relative to the other sites. Since a shift of several 
seconds is small compared to the time for drift in the av- 
erage trigger rate, the randomized data should reproduce 
all aspects of the real data associated with accidental co- 
incidences. Deviations of correlations observed in the real 
data relative to the randomized data could be indications 
of real correlations (i.e., not accidental) in the data. 

During the 17 month period, 1.68 x 10 7 trigger events 
were selected (according to criteria described below) and 
analyzed for time correlations. The distribution of time 
differences between consecutive triggers falls exponen- 
tially as expected, but does not follow an exact exponen- 
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FIG. 3: Probability of an excess trigger coincidence per site 
plotted vs. cumulative coincidence time scale. The upper 
plot shows the results for short time scales r < 10 -3 s and 
the lower plot shows the longer time scales 10 < r < 1 s. 
The lines drawn are the upper end (solid), center (dotted) and 
lower end (solid) of a 90% Feldman-Cousins confidence inter- 
val. The upper curve is interpreted as an upper limit for the 
probability per site of real coincidences for the corresponding 
time scale. 



to 
a. 

i— 

<u 




250 300 
At (US) 

FIG. 4: Trigger pair time distribution for At < 300 /is. Plot- 
ted points with error bars are real data and the horizontal dot- 
ted line indicates the expectation from random coincidences. 



tial distribution due to the varying number of operating 
detector sites. For a coincidence time interval At we 
count the number of consecutive trigger pairs with time 
separation less than At and subtract the corresponding 
number of pairs in the randomized data to form an excess 
-^excess) which may be positive or negative. If both mem- 
bers of a successive pair are from the same site, that pair 
is not counted (to eliminate instrumental effects such as 
PMT afterpulsing) . Fig. [21 shows N cxccss as a function 
of the coincidence time At. The large signal correspond- 
ing to real air showers is evident at short time intervals 
At < 2 (is. If 3 of the 60 sites are omitted (see Fig. 1), 
all pairs with separation 1 km or less are removed from 
the data set and the air shower signal disappears. In the 
following, we select sites so that relative distances are 
always greater than 1 km. 

In order to search for correlations on any time scale 
less than 1 second, we compute the probability of an ex- 
cess coincidence per site for the cumulative time interval 
{0, t} according to 
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in which (iV s it es ) = 37 is the average number of opera- 
tional sites and iVtrig is the total number of trigger pairs. 
We compute a 90% confidence interval for P CKC css(t) us- 
ing the method in [Tlj . and interpret the upper limit as 
the 90% confidence level upper limit for the excess prob- 
ability per site for the interval {0,t}. The results are 
shown in Fig. [3] 

The upper curve in Fig. should grow as the y/r, and 
it does increase in a fashion that is consistent with ex- 
pectations. One expects that the lower 90% C.L. should 
behave in a symmetric manner, and remain negative or 
close to zero. However, Fig. [3] indicates a rather substan- 
tial excess at time scales of order ~ 10~ 45 s, or ~ 30 /is. 
The excess heals itself at larger r as we add much more 
data that shows (apparently) no correlations. A less sig- 
nificant excess is perhaps evident at ~ 0.1 s. 

A histogram of time differences for real data is plotted 
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FIG. 5: Time and distance distributions for coincidences sep- 
arated by less than 50 /ts. Real data is shown as points with 
error bars, randomized data by dotted lines. The distance 
distributions for both the real and randomized data are re- 
lated to the spatial distribution of detector sites in Fig. 
but the difference between the real and random histograms 
potentially contains information on the excess coincidences. 



along with random expectations for short time scales in 
Fig. 01 which indeed shows an excess in the — 50 fis 
bin as one might expect from Fig. [3] The random data 
are averaged over 1 ms and plotted as a constant line to 
improve the statistics. The excess in the first bin is 2.5a 
over the expectation, significant at 99.26% C.L. 

In Fig. we display the data for the first bin of Fig. 21 
i.e., At < 50 /is. Although there is a preference for 
shorter time differences < 10 /is, no clustering at short 
distances is evident in the spatial distribution, indicating 
the excess events are likely not due to single air shower 
events. The excess in the first bin at At < 10 fis is 
significant at 99.70% C.L. 

The excess coincident trigger events (relative to ran- 
dom coincidence expectation) within the — 50 fis time 
bin appear to be distributed over distances larger than 
10 km. The average separation distance for events in 
the 0—50 /is time bin is Ldata = 23.1 ± 0.7 km, com- 
pared to -Lrandom = 22.6 ± 0.16 km for the randomized 
distribution. Thus the apparent excess events exhibit an 
average separation distance of L excess = 27 ± 7 km, con- 
sistent with i ran( jom- I n contrast, we note that single air 
showers generated by primary cosmic rays near vertical 



incidence are much smaller, extending only over several 
km, and should exhibit a very steeply falling distribution 
in separation distance. 

Since the apparent excess in our data is associated with 
average coincidence time At ~ 17 /is and L ~ 27 km 
these events would be consistent with rather vertical 
incidence 9 ~ 15° for the primaries (or pairs of pri- 
maries). Highly inclined or horizontally incident single- 
shower events would imply much longer average At = 
L/c ~ 90 us than indicated by the data in Fig. 4. 

Finally, we note that we also attempted to find cor- 
relations of single isolated "match" events with distant 
trigger events. Due to the high rate of accidental "match" 
hits in that analysis, the results were much less sensitive 
and so are not reported here. 



Summary and Conclusion 

A search for time correlations in cosmic ray data col- 
lected by the CHICOS project has been performed. The 
results are generally consistent with a lack of any real 
correlation between isolated events, except for trigger 
events (E > 10 14 eV) with coincidence times At < 50 fis 
where an excess number of coincidences with a signifi- 
cance of 2.5cr is observed. This excess is distributed uni- 
formly over the 17 month data period, and the events 
are randomly distributed over the array sites within the 
400 km 2 area of the CHICOS array. For smaller times 
At < 10 fis the data indicate an even more significant 
excess of 2.9<r. We know of no previous cosmic ray ex- 
periment that would be sensitive to correlated shower 
events at these energies separated by 10 — 50 km. The 
LAAS experiment Q has only 8 sites at much larger dis- 
tances (~ 500 km) and would have seen < 1 event given 
the rate we observe (assuming the correlations persist to 
those larger distance scales). 

If the apparent excess of correlated trigger events at 
At < 50 /is is interpreted as a real signal, then we ob- 
serve a coincidence probability of P e xcoss ~8x 10~ 8 per 
site for each trigger signal. This low probability implies 
that observing triple coincidences is extremely unlikely. 
However, normalizing for the aperture of the detector 
arrangement at each site, one would then infer that a 
substantial fraction, perhaps 10 -3 or more, of high en- 
ergy cosmic rays come in pairs or multiplets within a 
~ 1000 km 2 area. However, this interpretation is strongly 
energy dependent, and since we have not measured the 
event energies this could be true over any subset of ener- 
gies in the range 10 14 — 10 18 eV. 

Correlated events at ~ 27 km separations would prob- 
ably result from rather local events, on the distance scale 
of the solar system. Sources at galactic distance scales 
(or greater) might generate correlated events over much 
larger distance scales, but it is extremely unlikely that 
particles could be reliably propagated at such short trans- 
verse distances over ~ 10 17 km. The photodisintegration 
of heavy nuclei (i.e., iron) by solar photons is a possibility 
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that has been previously considered [3 , but the predicted 
fluxes are extremely low. Based on the estimates in 0, 
we expect that about 10~ 5 of the incident cosmic rays at 
£~6x 10 17 eV will arrive as pairs within the size of the 
CHICOS array due to photodisintegration. These pairs 
would typically consist of a heavy nucleus along with a 
~ 10 16 eV nucleon. The rate of these pairs within the 400 
km 2 area of CHICOS would be about 1 per year. The 
efficiency for the ~ 10 16 eV nucleon to generate a trigger 
at any of the 60 CHICOS sites is quite small, ~6x 10~ 3 , 
so the apparent excess coincidence rate (~ 40 events per 
year) is about 10 4 times greater than the expected rate 
for this process. Thus if wc maintain that the corre- 
lations originate in the solar system, we probably need 
either anomalous interactions with the solar wind or with 
solar photons. 

The distance distribution of excess events in Fig. [5] 
yields an average separation distance L exccss — 27±7 km, 
much larger than extensive air showers from single pri- 
mary particles near vertical incidence. Although it seems 
unlikely, it may be that the excess coincidences are asso- 
ciated with subtle unknown properties of large air shower 
events. It may be possible to explore such a possibility by 
generating unthinned simulations with high spatial reso- 
lution (~ 1 m) of air showers over a large ~ 1000 km 2 
area. Although this is a formidable computational task 
for present simulation programs such as AIRES [10|, we 
will attempt to address this issue in the near future. An- 
other unlikely possibility is that the generation of corre- 
lated shower events in the upper atmosphere separated 



by > 10 km at ground level could signal the onset of 
exotic new phenomena. 

It is clearly desirable to obtain more data to improve 
the statistical precision and attempt to verify or refute 
this observation. We hope to continue to expand and 
operate the CHICOS array to obtain a larger data set 
over the next few years. In addition, it may be possible 
for the Auger observatory to search for such correlated 
events. Further simulation studies of the properties of 
extensive air showers at large distances would also be 
helpful in the interpretation of these data. 
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